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Abstract 

Second generation high energy neutrino telescopes are being built to reach sensi- 
tivities of neutrino emission from galactic and extragalactic sources. Current neu- 
trino detectors are already able to set limits which are in the range of some emis- 
sion models. In particular, the Antarctic Muon And Neutrino Detection Array 
(AMANDA) has recently presented the so far most restrictive limit on diffuse neu- 
trino emission (Achterberg et al. (IceCube Coll.), 2007). Stacking limits which apply 
to AGN point source classes rather than to single point sources Achterberg et al. 
(IceCube Coll. &; P. L. Biermann) are given as well. In this paper, the two different 
types of limits will be used to draw conclusions about different emission models. An 
interpretation of stacking limits as diffuse limits to the emission from considered 
point source class is presented. The limits can for instance be used to constrain the 
predicted correlation of EGRET-detected diffuse emission and neutrino emission. 
Also, the correlation between X-ray and neutrino emission is constrained. Further 
results for source classes like TeV blazars and FR-II galaxies are presented. Start- 
ing from the source catalogs so- far examined for the stacking method, we discuss 
further potential catalogs and examine the possibilities of the second generation 
telescopes IceCube and KM3NeT by comparing catalogs with respect to north- 
ern and southern hemisphere total flux. 
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1 Introduction 



As is well known for hundreds of years, the sky would be infinitely bright 
in optical light, if all space were homogeneously full of stars , and space were 
infinitely extended without any change in properties (jOlbersl . Il826l ). This rea- 
soning can be applied to stars as well as to galaxies. This argument can also 
be applied to electromagnetic emissions and just as well to neutrino emission. 
The solution to Olbers paradox is that space is not homogeneously full of 
identical distributions of stars or galaxies, and that the Universe evolves quite 
strongly. This implies that the integral over the optical light of all galaxies has 
a finite sum, which does not exceed an observable level. Clearly, given a proper 
sampling of galaxies, it is possible to determine this level. In some cases, there 
is a unique relationship between the optical light and the emission at some 
other wavelength. Starting from this correlation, the summed emission from 
all galaxies over all history can be deduced and it is a simple function of the 
sum at optical wavelengths. The optical emission can certainly be replaced 
with the emission at some other electromagnetic wavelength, such as X-rays. 
Equally, the correlated wavelength can be neutrino emission. If we now do not 
know what the sum of the emission is at the electromagnetic wavelength, but 
have limits, then a limit for the neutrino emission can be deduced a fortiori. 
Vice versa, if we have a limit for the neutrino emission, a corresponding limit 
for the electromagnetic emission can be derived. In either case, it might hap- 
pen, that we have a known background, and then can deduce whether this 
specific class of sources could possibly explain all background at the other 
"wavelength". We start this argument with neutrinos, set a limit, and then 
derive a limit for the electromagnetic background. We could also start with 
the electromagnetic background, and then set a limit for the neutrino emis- 
sion. So we have an Olbers paradox for neutrinos. This is the key point of this 
paper. 



More concretely, the interpretation of different AMANDA neutrino flux lim- 
its is done for different classes of Active Galactic Nuclei (AGN). In the case 
of a stacking analysis in which the strongest sources from the same class are 
selected, a method for the interpretation of these limits as limits to the diffuse 
flux from the given source class is developed. Additionally, future possibilities 
of second generation neutrino telescopes like IceCube and KM3NeT are 
examined. Different source classes are discussed with respect to their distri- 
bution in the sky, i.e. which sources are in the northern and which are in the 
southern hemisphere. The field of view for IceCube respectively KM3NeT 
is the northern respectively the southern hemisphere and their view of the sky 
is complementary. 
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This paper is organized as follows: in Section 2, the prevailing models are 
discussed according to the normalization options given by experimental ob- 
servations of Cosmic Rays (CRs) and photons at different wavelengths. It is 
essential to discuss in this context which stacking limit can be applied to 
what prediction. Current neutrino flux limits are investigated in Section 3. 
A general ansatz for the interpretation of stacking limits as diffuse limits is 
presented here. In Section 4, the method is applied to AMANDA's limits 
which have been derived by stacking AGN according to their electromagnetic 
output. The question which stacking limits apply to which normalization sce- 
nario is discussed in Section 5 and it is examined how the limits constrain 
different models. In Section 6, further source classes are examined according 
to the possibility of applying the stacking method for high energy neutrinos, 
i.e. E v > TeV. In addition, a comparison of the contribution from the northern 
and southern hemisphere is done for various source classes in order to compare 
the capabilities of next generation's telescopes IceCube and KM3NeT. The 
results are summarized in Section 7. 



2 Discussion of prevailing neutrino flux models 



In some hadronic acceleration models it is assumed that for each class of AGN, 
the electromagnetic emission is correlated to a neutrino signal. Apart from in- 
dividual normalization factors, the corresponding cosmological integrations are 
basically mathematically identical. In this section, the correlation between the 
emission of neutrinos and photons at different wavelengths will be discussed 
according to neutrino flux models which are currently being discussed in the 
literature. 



2.1 TeV photon sources and Cosmic Rays 



Sources of electromagnetic TeV emission c an be interpreted as optically thin to 



photo n-neutron interactions, r in <C 1 (e.g. lMiicke et al.l . l2003l ; lMannheim et al. 



20011 ) in hadronic acceleration models. In such a scenario, charged Cosmic 



Rays (CRs) are produced in the vicinity of the source through the decay of 
the escaping neutrons. In this case, the resulting neutrino energy density would 
be proportional to the extragalactic CR component measured at Earth. A the- 
oretical upper bound of such a contribution to the diffuse neutrino flux has 
been derived by Mannheim, Protheroe and Rachen (in the following referred 
to as the MPR bound), see ( Mannheim et al. . 200 ll ). Within the framework 
of the proton-blazar model, the neutrino flux from High-frequency peaked BL 
Lacs (HBLs) has been calc ulated using the con nection between Cosmic Rays 
and neutrino emission, see (IMiicke et all 120031 ) . 
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2.2 Normalization to the diffuse EGRET and COMPTEL spectrum 



In the case of optically thick sources (r 7 „ > 1), the photons interact with nu- 
cleons in the source before escaping at lower energies leading to the emission 
of sub-TeV photons. Therefore, the di ffuse extragalacti c background deter- 
mine d by the EGRET0 experiment (ISreekumar et all Il998l ; IStrong et al.l 



20051 ) (i? 7 > 100 MeV) can be interpreted as an avalanched TeV signal from 
blazars and can thus be used to normalize the neutrino flux from EGRET 
or COMPTElUl-type sources. Again, MPR give an upp er bound to the con- 



tribution from such sources in (IMannheim et al.l . 120011 ). which is much less 



restrictive than the optically thin case. Apart from the bound, a calculatio n 
of the maximum contribution from blazars is given in (IMannheim et al.l . l200ll ). 
In addition to the contribution from HBLs within the framework of the proton- 
blazar model, a contribution from the optically thick Low-frequency peaked 
BL Lacs (LBLs) can be calculated using the EGRET diffuse extragalactic 



backg round for a normalization of the neutrino spectrum, see (IMiicke et al. 



20031 ). A model of p'j interactions in AGN and collis ions of protons fro m the 
core with protons of the host galaxy is derived in (jMannheiml . Il995l ). This 
model will be referred to as M95-A. Alternatively to optically thin sources 
with photon emission above 100 MeV, the environment can be optically thick 
allowing only photons below 100 MeV to escape. In that case, the neutrino 
signal can be normalized to the diffuse extragalactic cont ribution measured 



by C OMPTEL at energies in the range of (0.8, 30) MeV flKappadath et al. 



19961 ). This would enhance the contr ibution of neutrin os from p'j by almost 
an order of magnitude as shown in flMannheiml . Il995f ). We will refer to this 
model as M95-B in the following. 



A model by lStecker fc Salamonl fll996h was originally using the diffuse cosmic 
X-ray background (see Sec. 2.3) and has recently been modified in a way 
that it is u sing the COM PTEL diffuse background to normalize the neutrino 
spectrum ( Steckerl . 2005 ). This reduces the formerly very high contribution 
by a factor of 10. In addition, oscillations have been taken into account which 
leads to a further reduction by a factor of 2. 



2.3 The ROSAT X-ray background as normalization option 



The measurement of the diffuse extragalactic contribution in X-rays by RO- 
sat3 has raised the question whether it is produced by radio-weak AGN. As- 
suming that the X-ray emission comes from the foot of the jet, the X-ray signal 

1 Energetic Gamma Ray Experiment Telescope 

2 COMpton TELescope 

3 ROentgen SATellite 
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would be acco mpanied by a neutrino flux. A model by lNellen et al.l (119931 ) and 
an approach by lStecker &: Salamonl (119961 ) have been presented. An alternative 
scenario would be the up-scattering of thermal electrons via the Inverse Comp- 
ton effect. In that case, the X-ray component w ould not be accompanied by a 



neutrino signal, see e.g (IMannheim et al.l . ll995l ). U ntil today, about 75% of th e 



diffuse X-ray signal has been resolved by_ti 
with the help of Cha ndra and XMIVp- Newt on data 



ROSAT ( IBrandt fc Hasingerl . 120051 ). 

this number can be 



updated to 90%, see (IBrandt fc Hasingerl . 120051 ) and references therein. More 
than 70% of the diffuse background could be connected to the X-ray emis- 
sion of Active Galactic Nuclei most of which are radio weak. AMANDA's 
measurements of a diffuse neutrino flux did not yield a significant signal and 
constrains these models strongly. 



2-4 Correlated radio and neutrino emission 



Detailed examination of multi-wavelength observations have sho wn a correla- 



tion between the disk and j et power of Active Galactic Nuclei, see (IFalcke et al 



19951 ; IFalcke &: Biermannl . Il995l . Il999l ). This correlation can be used to deter- 



mine the diffuse neutrino flux from FR-II galaxies and flat s pectrum radio 



quasa rs within the framework of the jet-disk symbiosis model ( iBecker et al. 
20051 ). The neutrino flux has been assumed to follow an E~ 2 spectrum which 



leads to a relatively high neutrino flux at energies above TeV energies. As- 
suming on the other hand a correlation between the radio spectral index and 
the index of the proton spectrum leads to an E~ 2S neutrino spectrum. Since 
the flux is normalized at low energies, the contribution decreases significantly 
at higher energies with an increasing spectral index. 



3 Neutrino flux limits 



The conversion of stacking limits into stacking diffuse limits is the main topic 
of this section. Basically we use the contribution of the identified sources to the 
integrated background at a chosen electromagnetic wavelength to estimate by 
way of a physical model the corresponding ratio for neutrinos: What fraction of 
identified neutrino source candidates goes towards the integrated background? 
It is important to note here that we use conservative estimates in order to ob- 
tain absolute, reliable upper limits. As limits to the flux of identified source 
classes we use the sampl e of AGN classes from the AMANDA stacking analy- 



sis which is described in (|Achterberg et al. (IceCube Coll. &: R L. Biermann) 



20061 : lAchterberg et al. (IceCube Colli l2007f ). 



X-ray Multi-Mirror Mission 
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To avoid any confusion of the notation, a paragraph on the limit conventions 
precedes the actual discussion of the limits. 



3.1 Limit Conventions 



Throughout this paper, several types of neutrino flux limits will appear. Limits 
are usually given in the form E 2 ■ dN/dE with dN/dE as the differential 
neutrino flux at Earth and E as the neutrino energy. 

Throughout the paper, the limits will be denoted as follows: 

• <f>DL. I3iff use Limit (DL) derived by using data from the complete northern 
hemisphere. It is given in units of GeV cm _2 s~ 1 sr _1 . 

• Stacking Limit (SL) in units of GeV cm _2 s~ 1 , obtained for the point 
source flux from a certain class of AGN. The principle is indicated schemati- 
cally in Fig. 1. As an example, three source classes, FR-I and FR-II galaxies 
as well as blazars, are displayed. The stacking limits are obtained by stack- 
ing the most luminous sources of the same class in the sky, indicated by 
sources with filled circles. Weaker sources (empty circles) are not included. 

• Q D1J1J : Stacking Diffuse Limit (SDL), derived from the stacking limit in 
the same units as the diffuse limit, GeV cm _2 s~ 1 sr" 1 . It is determined by 
taking into account the contribution from weaker sources as well as yet 
unidentified sources, present in a diffuse background. 

A similar convention is used to denote the corresponding sensitivities: 

• $ UD : Diffuse Sensitivity (DS), giving the sensitivity of AMANDA towards 
a diffuse signal from the northern hemisphere. Units are GeV cm -2 s" 1 sr _1 . 

• $ : the sensitivity to a single point source in AMANDA. Units are 
GeV cm -2 s" 1 . 



3.2 Limits on the total diffuse neutrino flux 



The unfolded spectrum of the diffuse neutrino flux as observed by the 
AMANDA experiment is indicated in e.g. Fig. 8. It can be seen that the 
measured flux follows the prediction of the atmospheric neutrin o spectrum and 
no extragalactic contribution is observed (IMunich et al.l . 120051 ) . The unfolded 
spectrum is derived using data from the year 2000. The current limit on the 
diffuse muon neutri no flux from extragalactic sources isgiven for four years 
of data (2000-2004) flAchterberg et al QceCube Colli l2007h . 



DL 



7.4- 10- 8 GeVcm-V 



sr 
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Fig. 1. Schematic figure of the principle of source stacking. Objects are clas- 
sified into source classes (here: FR-I and FR-II galaxies as well as 
blazars). Under the signal hypothesis that photons in a certain band- 
pass are directly correlated to neutrinos, the sources are selected ac- 
cording to their strength in the corresponding bandpass. The most 
intense sources are then stacked for an analysis of a neutrino signal 
from the source class. This is indicated for the case of FR-I galaxies 
as circles around the objects. Only the strongest sources are used in 
the stacking analysis (filled circles), while weaker objects are not used 
(empty circles) . An estimate of the remaining source signal is needed 
as well as knowledge about unresolved sources, possibly seen as diffuse 
background in order to interpret the stacking neutrino limits as diffuse 
upper bounds to the given source class. The diffuse background varies 
with the source class. 



in the energy range of 10 4,2 GeV to 10 6 ' 4 GeV. 

In addition to a diffuse search, AMANDA examines a possible signal from 
single sources. No signal excess above the atmospheric background was ob- 
served yet. To maximize the significance of the signal to background ex- 
pectation, a stac king method has been applied examining different sourc e 



classes of AGN (lAchterberg et al. (IceCube Coll. fc P. L. Biermannl 120061 ) 



This method is commonly used in cases where a single source is not l ikely to 
have a significant signal above background (e.g. lMannheim et al.l . ll996l ). Given 
the limit for a certain source class, $ , the differential flux renormalized to 
a diffuse signal is 



SDL 



27rsr 



(2) 



The transition from point source to diffuse flux is done by dividing the source 
limits by AMANDA's field of view, 2 7r sr. The stacking factor e is the ratio of 
the total photon signal in the sample and the photon signal which is included 
in the analysis. Figure 1 gives a schematic representation. The stacking factor 
corrects for all sources with empty circles in the diffuse limit, since the signal 
from these sources is not included in the stacking analysis. The diffusive factor 
£ represents the ratio of the total diffuse flux to the contribution of resolved 
sources. Sources which have not been resolved yet have to be considered as 
potentially contributing to the total diffuse neutrino signal as they contribute 
to the isotropic photon signal. In Fig. 1, this is indicated by the diffuse back- 
ground which is different for each source class. This background is identified 
in a few cases by looking at the luminosity function of the objects, in others 
by measurements of the diffuse component. A detailed view on e and £ will be 
given in Section 4. 



4 Point source limits on AGN neutrino fluxes 



Recently, a stacking analysis has been published for the first time for a neutrino 
signal from Active Galactic Nuclei for which eleven AGN source classes have 
been defined. The source classes have been defined and selected as follows. The 



catal ogs below have been selected in (lAchterberg et al. (IceCube Coll. fc P. L. Biermann) 



2006|). 



(1) GeV blazars detected by EGRET (GeV). 

(2) Unidentified GeV sources observed by EGRET (unidGeV). Sources of 
possible galactic origin have been excluded. 

(3) Blazars that have been observed at infrared wavelengths {IR)- 

(4) Blazars that have been observed at keV wavelengths by the experiment 
HEAO-ALI] (keV(H)). 

(5) keV blazars observed by the ROSAT experiment (keV (R))- 

(6) The class of blazars with observed TeV emission (TeV). 

(7) Compact Stee p Sources (CSS) and GH z peaked sources (GPS) (CSS/GPS) 



as selected by lO'Dea &; Bauml (119971 ) . 



High Energy Astronomy Observatories-A 
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Fig. 2. Limits on the neutrino flux for a given source class. Limits in this 
representation are given per source. The source classes are labeled 
with the abbreviation indicated in the list of the samples above. 



(8) Low luminosity Faranoff Riley galaxies (FR-I) with M-87 included (FR- 
I(M)). 

(9) FR-I galaxies excluding M-87 {FR-I). These two different cases are nec- 
essary to consider, since M-87 is the closest AGN (distance ~ 20 Mpc) 
and dominates the total signal of all FR-I galaxies. 

(10) High luminosity Faranoff Riley galaxies {FR-II). 

(11) Quasi Stellar Objects (QSO). 



References to the corresponding catalogs as well as a detailed description of 
the analysis methods are given there. In brackets, the abbreviations for the 
source classes are given as they appear in Fig. 2. 



The limits on the neutrino flux from the given source classes are displayed 
in Fig. 2. They are compared to AMANDA's single source sensitivity of 
$PS = 5.9 . iQ-8 GeV cm s to point out the improvement that has been 
achieved in the stacking approach. In the following, the stacking factor e and 
the diffusive factor £ will be discussed. It will be shown that seven of the eleven 
samples can effectively be used to determine a diffuse limit on the neutrino 
flux from the given class. 
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4-1 The Stacking factor 



The number of sources in a sample that is included in the stacking analysis is 
determined by optimizing the expected significance for signal and background. 
Subsequently, the amount of the total photon signal of the analyzed sources is 
determined by the number of included objects. This is shown schematically in 
Fig. 1. The strongest FR-I galaxies (filled circles) are included in the stacking 
analysis. The stacking process for FR-I galaxies is indicated by circles around 
the stacked objects. FR-I galaxies having a too weak signal (empty circles) are 
not part of the stacking and have thus to be included in the stacking factor 
e. Quantitatively, the ratio of the total photon signal and the signal of the 
sources that are used in the stacking approach is represented by e. The precise 
values of e for each class are given in Table 1. 

Some of the stacking classes cannot be used for a diffuse interpretation. Firstly, 
the number of sources needs to be sufficiently high. An indication for an in- 
sufficient number of sources would be the use of all the sources in a sample 
in the stacking analysis. Secondly, the luminosity function should be moder- 
ately falling: if only very few sources in our vicinity dominate the signal, no 
conclusions about the overall diffuse contribution can be justified. 

The complete sample was used for stacking in several cases due to the low 
number of detected sources in the given wavelength interval. In particular, 
these source classes are the infrared and HEAO-A sources as well as the TeV 
blazars. Since all twelve observed infrared sources have been used and only 
three sources are reported from HEAO-A, a diffuse limit to the source class 
cannot be derived for these samples. There are five sources in the stacking 
sample of TeV blazars. Due to the small number, this class cannot be used 
for a diffuse interpretation. Due to the strong absorption of the TeV signal in 
photons, the contribution of a neutrino flux from TeV-resolved photon sources 
is expected to be much lower than the total flux. This is discussed in detail in 
section 5. 



Apart from the three classes already excluded, the class of FR-I galaxies in- 
cluding M-87 is not suitable for an interpretation of a stacking limit as a diffuse 
limit: M-87 makes up most of the photon signal from all FR-I galaxies. This is 
also the reason why the stacking limit is very close to the single point source 
sensitivity. This result must therefore rather be interpreted as a point source 
result of M-87. The FR-I sample without M-87 on the other hand is applicable 
for the diffuse interpretation. 
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no 
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178 MHz 


1.1 
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2.91 
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178 MHz 


2.65 
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30.4 


yes 


QSOs 


UV 


1.3 


model dep. 


6.70 
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Table 1 

Table of the source class limits obtained with the stacking method. Five 
years of data, 2000-2004, have been u sed for the analysis with AMANDA in 



(jAchterberg et al. (IceCube Coll.)l . 120071 ). Listed are the source class, the selection 
wavelength, the stacking and diffusive factors, the stacking limit <£ OLj in units of 
10~ 8 GeV cm -2 s" 1 and the possibility of interpreting the stacking result as a diffuse 
limit. 

4-2 The diffusive factor 



Apart from the contribution of resolved sources to a diffuse background, 
a component of unresolved sources has to be considered. The ratio of the 
total diffuse signal to the signal of resolved sources is called the diffusive 
factor £ in the following. In Fig. 1, this is indicated as the diffuse pho- 
ton background of unidentified sources. A truly diffuse contribution from 
e.g. dark m atter decays is discussed as a component of the diffuse photon 
background (IBiermann fc Kusenkd . 120061 ). but it will be considered as negligi- 
ble here. Any contribution of a diffuse component not connected to any source 
population would improve the limits accomplished. 



70% of the diffuse X -ray background has bee n identified as AGN in the 
ROSAT data sample (IBrandt fc Hasingerl . 120051 1 . Thus, the diffusive factor £, 
which is the inverse of the fraction of resolved sources, for the keV ROSAT 
data sample is given as £rosat = 1-43. Note that today, XMM-Newton 
and Chandra provide values of ~ 95% resolved sources. But since both the 
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neutrino flux predictions are based on and the stacking limit has been derived 
from ROSAT data, a fraction of 70% resolved sources is used. 



For EGRET sources, it could be shown that only 1/12 of the total diffuse 
backg round can be ma de up by resolved sources, as d erived by lChiang &: Mukherjee 
fll998f ). The claim from lChiang & Mukherieel fll998h that about 25% of the dif- 
fuse background of blazars is made up by resolved AGN, while the remaining 
contribution is from unresolved sources which are not blazars, wo uld reduce 
£ significantly to Regret ~ 4. However, IStecker &: Salamonl (120011 ) point out 
that this estimate is based on assumptions that do not hold for EGRET 
blazars. Thus, £ max = 12 will used as a conservative number in the following. 
With the launch of glastOJ the question of the contribution from resolved 
blazars will most likely be settled. 



For the remaining source classes, there is no explicit estimate of the contribu- 
tion of resolved sources to a diffuse background. In this case, an estimate of £ 
can be possible for source populations where there is an estimate of the weakest 
sources in the sky. In this case, the luminosity evolution of the population has 
to be considered. In the Euclidean case, sources evolve as N(> S) oc S~ 3 ^ 2 with 
S as the observed flux and N(> S) as the number of sources for fluxes greater 
than S. In such a case, the integral of the distribution diverges for a vanishing 
minimal source strength, N(S min — > 0) — > 0. For cosmological sources, the be- 
havior is not Euclidean an d in some case s, it is possible to determine a lower 
threshold for the flux, see ( Longairl . 1998 . e.g.). For instance, such determina- 
tion can be done in the case of FR-II galaxies which are high luminous AGN 
with a defined luminosity limit of L ra( ao ~ 10 43 erg/s (jFaranoff fc Rileyi . Il974l ). 



The determination of a lower luminosity limit enables the calculation of the 
total number of sources in a population N to t- The total number of sources is 
given by integrating the number per flux interval, dN/dS over the flux: 



tot 



dN 
~dS 



dS. 



(3) 



For the case of divergence for S m - m — > 0, the lower integration limit is essential 
to know as mentioned before. This is not trivial, since measurements for most 
source classes do not reveal the low luminosity cutoff. The luminosity functions 
of different AGN classes at z = are shown in Fig. 3. 

In this case, the diffusive factor can be estimated by assuming that sources 
which are not in the catalog are not stronger than the weakest source in the 
catalog. Thus, the flux of any additional source S l add is smaller than the flux 



Gamma-ray Large Area Space Telescope 
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Fig. 3. Radio Luminosity Functions at z = for Flat Spec trum Radio 
Quasars (FSRQs, solid line) ( Dunlop fe Peacock . 199C)h. FR-I and 



FR-II galaxies (dashed and dotted lines) dWillott et all bOQlh . 



as 

well as for QSOs (dot-dashed line) ISchmidtl (|l972l ). Density units 
are 1/Gpc 3 • AlogL and luminosity units are L42 := L/(10 42 erg/s). 
Only FR-II galaxies show a cutoff at low luminosities. All other source 
classes show an increasing or constant population with decreasing lu- 
minosity. 



of the weakest source, S wea k, S l add < S wea k for all non- resolved sources i. The 
total flux not considered in the calculation is thus 

Stot = J S dS < S wea k ■ N add . (4) 



Here, N aa -d = N to t — N cat is the number of additional sources which is calculated 
from the total number of sources expected in the sky, N to t, subtracting the 
number of sources in the catalog, N cat . Then, £ is calculated to 



qcat 1 qadd 
°to t ~r ° tot 
qcat 



£ tot 1 tot /r\ 

^~ qcat \°) 
'tot 



^tot "I" S W eak ' N a 
qcat 



(6) 



'tot 



where is the total flux in the catalog. 
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Figure 4 shows the behavior of the maximum value of the diffusive factor 
£max according to Equ. (6) at the example of EGRET sources. The figure 
emphasizes the challenge in such a representation. It is essential to have an 
estimate of the total number of sources in the class to get a good estimate 
for £ max . This implies the necessity of the knowledge of the absolute lower 
luminosity limit. For FR-II galaxies, this i s relatively simple, si nce the number 
of sources decreases for low luminosities ( Willott et all 1200 ll ). FR-I galaxies 
on the other hand have a high number of low luminosity sources and it is not 
known at which luminosity the function turns. 



= 



10 J rr 



10 Z r 
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Fig. 4. Maximum diffusive factor £ max versus total number of > 100 MeV- 
emitting sources for EGRET sources. With £ max = 12, the lower limit 
of the total number of sources contributing to the diffuse EGRET 
background is N to t > 3.9 • 10 3 . 



In such a case, the diffusive factor can be quite large, since a high contribution 
from low luminosity sources is expected. Thus, £ has explicitly been deter- 
mined according to the neutrino flux model which is tested by t he limits. For 
FR-II galaxies, the investigated neutrin o flux model p r esente d by lBecker et al. 
f l2005h uses the luminosity function bv IWillott et~al1 fl200ll h labeled "FR-II" 
in Fig. 3 . In this total number of ~ 10 sources is expected and the 

maximum diffusive factor is determined to £ max = 160. For the remaining 
source classes, FR-I galaxies, GPS/CSS and QSOs, the value of £ max will not 
be given, since no explicit neutrino model for these source classes is examined 
here. Note that £ max gives an absolute upper limit: every source not included 
in the stacking analysis is assumed to have the flux of the weakest included 
source, S wea k- This results in an overestimate of £, since it is likely that only 
a small fraction of the sources has such a high flux. 
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For EGRET sources, Equ. (6) can be used to determine a lower limit of 
sources contributing to the diffuse EGRET background. Figure 4 shows the 
dependence of the maximum value of £ on the total number of sources in the 
source class of blazars emitting at > 100 MeV. Using a diffusive factor of 
£ = 12, the total number of sources contributing to the total diffuse EGRET 
background is Nf f RET > 3.9 • 10 3 . 



4-3 Comparing stacking diffuse limits with overall diffuse results 
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Fig. 5. Diffuse 
solid 



limits of seven different AGN classes. The 
line represents the total diffuse limit of 



AMANDA (jAchterberg et al. (IceCube Coll.l 12003). Several of 



the stacking diffuse limits are more restrictive than the general diffuse 
limit. The dashed line indicates the Waxman/Bahcall (WB) bound 
which gives a theoretical estimation of the maximum contribution of 
AGN to the diffuse neutrino flux. 



Figure 5 shows the diffuse limits that could be derived from the stacking ap- 
proach for each source class. It is compared to the diffuse limit of AMANDA 
(solid line) for the operation time considered (2000-2003), 
$ DL = 7.4-10" 8 GeV cm" 2 s" 1 sr -1 , given in (jAchterberg et al. (IceCube Colh) 



20071 ). Several stacking diffuse limits are more restrictive than the general dif- 



fuse limits. The exact numbers are given in Table 2, where the stacking limit, 
the diffuse bound derived from the stacking method and the ratio of the latter 
and the general diffuse limit are compared. Note that in the case of CSS/GPS, 
QSOs and FR-I galaxies, the diffusive factor still has to be applied. 
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Figure 5 also indicates the Waxman/Bahcall bound (WB bound, dashed line). 
This bound gives a theoretical estimate of the maxi mum contribution to 



be ex pected in neutrinos from Active Galactic Nuclei (IWaxman &: Bahcall 



19991 ). This bound is valid for optically thin sources at an energy of E v 



10 eV, where it is most r estrictive. The description for all energies is given 



by iMannheim et al.l (120011 ). Five of the seven source class limits lie below 
the indicated WB bound. This shows that the sensitivity reached with this 
method in AMANDA is already extremely high. Next generation's neutrino 
telescopes such as IceCube and KM3NeT which will have a cubic kilometer 
of ice resp. water instrumented will therefore be able to provide more infor- 
mation either by giving very strong restrictions on prevailing models or by 
confirming a positive neutrino signal. 

Note that the class of TeV blazars is on the other hand an optimal candidate 
for a general diffuse search. The contribution of photon-observable TeV blazars 
is very small, since many sources are hidden in photons due to the strong 
absorption, while they would still be visible in neutrinos. This will be discussed 
in Section 5. 
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(Mannheim 1995) 
CMiicke et al.. 2003) 


unid. EGRET 


31.7 


6.56 


1.23 




infrared 


10.6 


EXCL.: 


too few 




keV (HEAO-A) 


3.55 


EXCL.: 


too few 




keV CROSAT) 


Q 71 


2 65 


2 79 


f KtprKP'r A7 S\a 1 a in nn 1 QQfi i 

fNellen et al 1993) 










( Alvarez-Mnniz Rr Meszaros 2004) 


TeV blazars 


5.53 


EXCL.: 


too few* 


(Mannheim et al.. 2001). 










CMiicke et al.. 2003). 










(Mannheim. 1995) 


CSS/GPS 


5.94 


1-23 • (model 


6.02/£ mO(feZ 




FR-I incl. M-87 


4.11 


EXCL.: M87 


dominant 




FR-I excl. M-87 


2.91 


0-51 • (model 


t^S/ (model 




FR-II 


30.4 


2.05 • 10 3 
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( Becker et al.. 2005) 


QSOs 


6.70 


1.39 • (model 
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Table 2 

Table of the source class limits obtained with the stacking method. Five 
year s of data, 2000-2004, have bee n used for the analysis with AMANDA 



in (jAchterberg et al. (IceCube Coll. )l . 120071 ). The stacking limit is given in 



units of 10 8 GeV cm 2 s 1 while the stacking diffuse limits is given in units 
of 10~ 8 GeV cm" 2 s" 1 sr _1 . * While the 

class of TeV blazars cannot be used to determine a stacking diffuse limit, the general 
diffuse limit gives an upper limit to the contribution of TeV-observable blazars to 
the total diffuse flux as it is shown in Section 5. 
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5 Direct implications for AGN neutrino flux models 



The diffuse limits discussed in Sections 1 and 4 constrain some of the currently 
discussed neutrino flux models. However, it is known that these models bear 
different uncertainties in both normalization and spectral shape due to a lack 
of knowledge of the conditions at the source. 

In this section, the neutrino flux models having been discussed in Section 1 
will be examined with respect to the limits obtained. 



5.1 TeV blazars 



The detection of TeV photon emission from distant sources is limited due 
to the ab sorption of high energy photons by the extragalactic background 



light (e.g. IStecker et al.l . Il992l ; iKneiske et al.1 . 120041 ) . The absorption factor r] 



describes the ratio of the total emitted TeV photon flux from HBLs and the 
TeV flux from HBLs up to a redshift z mSuX . It is a measure of the absorption 
of TeV photons. The ratio of a diffuse neutrino signal from photon-observable 
TeV blazars will be calculated using the general diffuse limit and taking into 
account the absorption factor for TeV photons: Only a relatively small fraction 
of all TeV blazars can be identified in TeV photons, since sources at high 
redshifts are hidden due to the strong absorption. 



In this paragraph, the absorption factor and the general diffuse neutrino 
limit are used to derive the maximum contribution of TeV photon-observable 
sources to the total diffuse neutrino flux. 



5.1.1 The absorption factor rj 



In the case of TeV blazars, the absorption factor t] is much greater than 
unity due to the strong absorption of TeV photons by the extragalactic back- 
ground light. HBLs seem to have no or eyen a slightly negative evolu t ion ( e.g. 



Beckmann et al.l . 120031 ; iBade et al.1 . Il998l ; lLaurent-Muehleisen et all Il999l ) . 



For a no-evolution scenario as it is discussed for BL Lacs, the co-moving 
density p(z) is considered to be constant with redshift, p(z) = constant. Using 
a negative evolution with less than (1 + z) -a2 does not change the results 
significantly. We also neglect the positive source evolution which has to be 
present up to a certain redshift z*: both effects positive evolution up to z* and 
negative evolution at higher redshifts are believed to cancel. Each effect for 
itself alters the result less than 10% in opposite directions. 
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The normalization of the co-moving density is not important in this calcula- 
tion, since only ratios of concrete values of p are considered and the constant 
of proportionality cancels. 

The neutrino flux from sources up to a certain redshift z max is given by 

bo- J E\z)- 2 .p(z)^dz, (7) 



dE 

2 = 



where E'(z) = (1 + z) • E is the energy of the neutrino at the source. The 
^-weighted flux is thus given as 

dE = 00 ' C(2max) • ( } 



Here, C( 2 max) = //=o x ( z + 1)~ 2 ' p( z ) dV c /dzdz is the evolution factor depending 
on the upper redshift integration limit z max . If all present sources are consid- 
ered, 2 max ~ 7 is a reasonable value. Recent searches for luminous galaxies at 
redshifts between z = 6 — 8 imply t hat only few ultra- luminous o bjects only 



exist beyond z = 7 (jive et al.l . 120061 ; iBouwens &: Illingworthl . 120061 ). At higher 



redshifts, no large galaxies are observed, leading to the conclusion that only 
smaller galaxies can be present. Also, the results to not change significantly 
when going up to higher redshifts: the major contribution comes from redshifts 
of z < 3. 

The absorption factor 77 is given by the ratio of the total flux and the contri- 
bution that is observed by TeV-photon experiments: 

^max) = 7 f^T • (9) 



Here, z max is the maximum redshift at which TeV photon sources can be 
identified by present Imaging Air Cherenkov Telescopes (IACTs). TeV sources 
have up to today been observed up to a redshift of z=0.186 - 1ES1101-232 
flAharonianl . l2006j l. As a conservative upper limit, the maximum redshift of 
observation with IACTs is taken to be z max = 0.3, since an upper bound is 
derived from these values. The absorption factor depending on z max is shown 
in Fig. 6. At z max = 0.3, the numerical value is 

r/(2 m ax = 0.3) = 54. (10) 



It is important to note that for a lower energy threshold (e.g. 30 GeV), IACTs 
could be able to detect sources up to z max ~ 1. This would lead to a much 
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Fig. 6. Absorption factor tj versus maximum redshift ^ m ax- 

For TeV 

sources, the most distan t source observed is 1ES1101-232 at 
z=0.186 (jAharonianl . l2006h . n{z max = 0.3) = 54 is used for TeV 
sources. 



lower absorption factor of 



??0max = 1) = 4.6 



5.1.2 An upper limit to the contribution from photon-observable TeV blazars 



Using the above reflections, one can argue that a diffuse limit can be used 
to estimate the maximum diffuse signal from HBLs resolved in TeV photons: 
From the existing general diffuse limit, it is known that no source type and, 
thus particularly also not the class of TeV blazars, can contribute more than 



DL 



7.4 



f (T 8 GeV cm 2 s' 1 sr~ 



E 



dN 
dE 



< $ 



DL 
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TeV 



An upper limit to the normalization constant 0o of the neutrino flux can be 
derived by inserting Equ. (8) with z max = 10 in Equ. (12): 



< 



cm 



(13) 
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The contribution of sources at z < z max is given as 



dN 
~d~E 



COmax) $ DL 



(^max) — 00 " C(^max) < 
TeV S,l iU J 



(14) 



with the final result of 
dN 



E 



dE 



1.$ 



DL 



(15) 



TeV 



The upper limit on the neutrino flux up to z max is shown in Fig. 7. The curve 
shows the no-evolution scenario. 




Fig. 7. Upper limit on the contribution of photon-observable TeV blazars for 
the no-evolution scenario for HBLs. 



The contribution of photon-observable TeV blazars to the diffuse neutrino flux 
is therefore limited to 



dN 
dE 



0.3) < 1.37- 10 J GeV cm^s 



sr 



(16) 



TeV 



where f](z mSuX = 0.3) = 54 is used as a lower limit. The contribution of sources 
observed by IACTs is thus about three orders of magnitude lower than the 
possible total contribution. This is displayed in Fig. 8 together with the general 
diffuse limit. The limit to the maximum contribution from photon-observable 
TeV blazars is shown (obs. TeV) as well. The contribution is a factor of r]^ 1 ~ 
0.019 lower than the diffuse limit. This indicates that a diffuse analysis of 
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TeV sources is most effective with an overall diffuse approach, since most of 
the TeV sources are hidden due to the strong absorption at such high photon 
energies. The stacking of the strongest sources in the sky in comparison only 
comprise a small fraction of the total diffuse flux as these calculations show. 
This result stands in contrast to the other sources samples where the selection 
of the strongest sources yields stronger restrictions than a general search for 
a diffuse signal, see following paragraphs. 
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Fig. 8. Neutrino predictions for optically thin sources. TeV photons from de- 
caying pions escape together with neutrons which make up charged 
CRs by decaying in the vicinity of the source in protons. All models 
have been c orrected for oscillations . The MPR-bound derivation is de- 
scribed in ( Mannheim et all l200lh . Predictions for the contribution 



from HBL sources (HBL-I and HBL-II) depend on the extension of 
the acceleration region and the jet frame target photon density. Model 



1 and 2 represent two different parameter settings, see (IMiicke et al 



2003|). The general diffuse limit is shown as well as the limit to photon- 
observable TeV blazars {obs. TeV), indicating that the major contri- 
bution to the diffuse neutrino signal from TeV blazars is expected 
to come from hidden TeV photon sources for two different maximum 
redshifts, z = 0.3 and z = 1. 



The importance of reaching down to low energies with IACTs is demonstrated 
in this calculation, since lowering the energy threshold to observe sources up 
to z — 1 would already include sources with a one order of magnitude higher 
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maximum neutrino flux of 
dN 



E 2 



dE 



Omax) < 1.60 ■ 1(T 8 GeV cm" 2 s" 1 sr" 1 . (17) 

TeV 



5.2 Optically thick cases: sources of MeV and GeV 7 emission 



If the sources are optically thick to photon-neutron interactions, r 7n >> 1, 
the photon signal which is emitted from the sources lies in the MeV to GeV 
range. Thus, EGRET and COMPTEL diffuse cosmic photon fluxes are used 
to determine the expected neutrino contribution. 



5.2.1 EGRET blazars 

Figure 9 shows neutrino flux predictions in connection with EGRET data 
and the derived stacking diffuse limits for a neutrino signal from EGRET 
sources. Since the indicated models are all normalized to the EGRET diffuse 
spectrum, the limit of identified EGRET blazars applies. The gray line be- 
tween 1 TeV and 10 2 ' 8 TeV indicates the most conservative calculation where 
it is assumed that the total diffuse background as measured by EGRET it 
produced by AGN (£ = 12). The limit falls short of the MPR bound and also 
touches the maximum prediction for the source class predicted by MPR and 
it is possible to reach sensitivities at the level of predictions concerning a neu- 
trino signal from blazars. The same is valid for the prediction by Mannheim 
(M(95), A): The proton-proton contribution from optically thick blazars can 
be constrained even though the atmospheric flux exceeds the prediction at 
these energies. With an overall diffuse analysis, it would not be possible to 
extract any information on this low-energy part of the spectrum from AGN. 

The contribution from LBL within the proton-blazar model dominates the 
total observable flux at ultra high energies. With the current data, it is not 
possible to get to any reliable conclusion from the stacking method, since it is 
valid at energies E < PeV while the flux is present at much higher energies. 
A stacking analysis of cascade events with second generation telescopes such 
as IceCube would be an option to explore regions of energies at > PeV. 
A stacking approach for cascades is however accompanied by the challenge of 
getting a reasonable directional reconstruction of the neutrino events, since the 
cascade signal is generally not as boosted as a signal from neutrino-induced 
muons. 

Even in the most conservative case of a diffusive factor of £ max = 12, the 
model of proton-proton energies can be constrained at low energies. With the 
stacking method, the diffuse limit can be extended to sensitivities far below 
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Neutrino flux predictions for optically thick sources, normalized to 
the diffuse EGRET flux above 100 MeV. Model 1 (MPR) repre- 
sents the maximum contribution from blazars within the framework 



of lMannheim et al.l (|200ll ) in which the MPR bound for optically thick 
sources (dashed line) is derived. Model 2 (M95,A) - gives the neu- 
tri no spectrum from pp and pj interactions in blazars determined 
Mannheim! (|l995l ). Models 3 and 4 (LBL-I and LBL-II) represent 



in 



the predictio n of neutrino s from LBLs given two different parameter 
settings, see iMiicke et al.l ()2003l ). The stacking diffuse limit for the 



12 



class of identified EGRET sources is indicated with £ — £ma: 
where it is assumed that AGN produce the complete diffuse EGRET 
signal. 



the flux of atmospheric neutrinos. 



5.2.2 COMPTEL blazars 



A stacking analysis using COMPTEL blazars has not been done yet. The op- 
tion is discussed in Section 6. Current diffuse neutrino flux limits are still about 
a factor of 2 or more above the predictions. The hypothesis of an avalanched 
TeV signal down to MeV energies will be tested by IceCube. 
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5.3 The diffuse X-ray background and Radio Weak AGN 



The X-ray component of radio weak AGN can be correlated to the emission 

of neutrinos at the foot of the jet. This has bee n investig ated in calculations 

oflNellen et al.l (119931 ). IStecker &: Salamonl (119961 ) as well as I Alvarez- Muniz &: Meszaros 
( 120041 ) . All models are constrained by the current general diffuse AMANDA 
limit. The resolved ROSAT sources which have been used to determine the 
stacking limit are radio loud and thus, the limit is not as easily applicable to 
these predictions. If the same production mechanism is assumed in radio loud 
as in radio weak sources, the stacking diffuse limit of ROSAT sources would 
apply in this case. 

CRT 

However, the limit to radio weak objects $ r q is about a factor of 10 
higher than the calculated limit for radio loud sources $ r u since radio 
loud objects are about a factor of 10 less frequent than radio weak sources, 



$SDL < 



1.72 ■ 10~ 7 GeV cm _2 s -1 sr- 



The more restrictive limit in this case is therefore the general diffuse limit. 



Figure 10 shows the prediction of neutrinos being produced coincident ally 
with X-rays at the foot of AGN jet s. The diffuse RQ SAT-measured back- 
groun d has been used in two models (INellen et all Il993l ; IStecker fc Salamon . 
1996 ) to normalize the ne u trino spectrum. In the third case presented by 
Alvarez- Muniz Sz Meszarosl (120041 ). the luminosity evolution function of radio 
quiet AGN has been used so that the result also applies to the correlation 
with the X-ray diffuse background. 



In the case of the prediction by lNellen et al.l (119931 ). AMANDA's diffuse neu- 
trino flux limit lies an order of magnitude below the flux. Stecker et al. updated 
their c alculation by n ormalizing to the diffuse flux as measured by COMP- 
TEL (jSteckerl . 120051 ). This is discussed in the context of the MeV photon 
emission in correlation with neutrinos in Section 5.2. The absence of neutri- 
nos from RO SAT-detected sources implies strongly that the X-ray emission 
from AGN cannot be directly correlated to neutrino emission. While the miss- 
ing emission of neutrinos disfavors a hadronic scenario for the X-ray emission 
from AGN, it is consistent with Inverse Compton (IC) models for the X-ray 
emission and underlines leptonic models for the production of X-rays in AGN. 



Calculations by Stecker fc Salamon ( 1996[) show a flux pea king at higher en- 
ergies c ompared to t h e pred iction from INellen et al.l (119931 ) . The main reason 



is that INellen et al.l (119931 ) use a more conservative estimate of the maxi- 
mum energy and a simpler approach for th e spectral behavior. The spectrum 
of the model by IStecker &; Salamonl (119961 ) is not exactly i?~ 2 -shaped, while 
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Fig. 10. IStecker fc Salamonl (jl99d ) (StSa96, model 2) predict a large contri- 
bution of neutrinos from radio weak s o urces using ROSAT data to 
normalize the spectrum. Nellen et al. ( 19931 ) (NMB, model 1) use 
the same normalization method. M odel 3 and 4 (AM-I and AM-II, 
(jAlvarez-Muhiz &: Meszarosl . 12004 )) use the luminosity function of 
radio quiet AGN and normalize the single source spectrum. While 
the black hole mass is varied in model 3, model 4 varies the accretion 
rate. All these models have been corrected for neutrino oscillations. 
All models are excluded in their published form by AMANDA's dif- 
fuse limit. The typical E~ 2 limit is represented by the blue horizon- 
tal line. The curved blue line, crossing the E~ 2 -limit represents the 
AMANDA limi t for the years 2000-2003 ca lcula ted for a flux s haped 
as the model by Stecker &: Salamon ( 19961 ) and Stecker ( 20051 ). 



the AMANDA limits are typically derived as suming an E 2 spectrum . The 
modeling of different spe ctral shapes, includin g IStecker fc Salamonl (119961 ). has 



been done recently, see (IHodges et all 120061 ) . The curved blue line ind icates 



the limit modeled according to the model of IStecker &: Salamonl (119961 ). The 
limit lies half an order of magnitude below the prediction. 

While the previously describe d models normalize the diffus e spec trum by using 
ROSAT data, calculations by lAlvarez-Muniz fc Meszarosl (120041 ) use the radio 
luminosity function for radio quiet AGN and a single source normalization. 
Further differences with respect to model 1 and 2 are the maximum energy 
and the break in the neutrino spectral shape. The underlying idea is the same 
by assuming proton-proton and proton-photon interactions at the foot of AGN 
jets, leading to X-rays in coincidence with neutrinos. This model overproduces 
neutrinos as well as the previous ones. 
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Three independent calculations on the correlation of X-ray emission from ra- 
dio quiet AGN and neutrino emission have been examined. The absence of 
neutrino emission from the foot of relativistic jets implies that the particles 
there are not accelerated to high energy, and then interact. However, the con- 
cept that the innermost ring of the accretion disk, just underneath the jet, 



temperature (e.g. 


Falcke et al.. 


1995: 


Biermann et al.. 


1995 




Rachen et al.. 


1995; 


Donea Sz Biermann. 


1996; 


Mahadevan . 


1998). is consistent with the re- 



suits. This implies that the spin parameter of the black hole is high, above 
0.95, and that hadronic interaction at weakly relativistic temperatures pro- 
duce charged pions, which decay and provide an energetic particle seed pop- 
ula tion for further acceleratio n downstream for the radio emission of the 
jet ( iGopal-Krishna et al.l . l2004j ). This is then consistent w i th an Inverse Comp- 



ton explanation of the X-ray emission (IMannheim et all Il995l ). and predicts 



that there should be a large production of energetic neutrinos at an energy 
commensurate with the pion mass. 



5.4 Radio galaxies 



The radio emission of Active Galactic Nuclei is likely to be d irectly corre- 
l ated to neutrino emission in the jet. This has been discussed by lBecker et al. 
( 120051 ) at the example of FR-II radio galaxies and flat spectrum radio quasars 
(FSRQs). In both cases, the normalization of the neutrino spectrum depends 
on many different intrinsic factors, i.e. the correlation between radio and disk 
luminosity and the parameterization of the maximum proton acc eleration en- 
ergy. T he jet-disk correlation of AGN has been worked out by iFalcke et al 
( 119951 ); iFalcke fc Biermann! (119951 . Il999l ). Although the model includes several 
parameters, all are fixed by the comparison with the data except the accretion 
rate. Within that model, it is shown that the parameter reach extreme values 
and are not strongly scattered. The jet-disk symbiosis model has been proven 
for different source types, reaching from microquasars to quasars. 



A parameter in the calculatio n of the neutr i no flu x is the optical depth of the 
source which is unknown. In (IBecker et all 120051 ) . the effective optical depth 



of the source, r e // is defined as product of the proton-photon optical depth 
r P7 and a reduction factor due to contributions from Bethe-Heitler production 



rf 



T, 



eff 



'PI 



It can be derived from the solution of the transport equation, the neutrino flux 
is proportional to r e ///(l — exp(— r e //)). For large depths, the neutrino flux is 
therefor proportional to r e jf. The effective optical depth has been chosen to 
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be T e ff — 1 in the calculations described above. The optical depth of proton- 
photon interactions is determined by the product of the photo-hadronic cross 
section a vl and the photon density in the source, n 7 . These in turn depend 
on the disk luminosity L disk and the extent r of the source such as 

r P7 = 800— , (19) 



with L46 := Ldisfc/lO 6 erg/s and ryr := r/10 17 cm. This shows that the photo- 
hadronic optical depth itself uncertain, since low values around unity are pos- 
sible as well as extremely high numbers. AMANDA limits will be used to 
derive limits on the optical depth of FR-II galaxies and FSRQs. 

Note that the determination of r e // happens only within the specific model of 
neutrino production as it is described above. The difficulty of drawing more 
general conclusions lies in the uncertainty of the spectral index. This highly 
depends on the original spectral behavior of the protons. This can only be 
simulated and is not directly observable. 



5.4-1 FR-II radio galaxies 

Using a generic E~ 2 spectrum for the neutrino flux from FR-II galaxies, the 
flux normalization has been derived to be 

$ = 1.43 • 10~ 7 GeV cnT 2 s^ 1 sr _1 . (20) 



Comparing the normalization in Equ. (20) with the diffuse limit shows that 
a flux with the chosen parameter settings is not detected. An upper limit to 
the optical depth can be derived from the limit yielding 

re// < 0.5. (21) 



Within three years of IceCube, it will be possible to explore sources with 
T e ff > 0.029. 

Figure 11 shows the the prediction of a neut rino flux from FR-I I radio galax- 
ies. Model (1) is the original calculation of iBecker et al.l (120051 ) modified by 
a factor r e jy=0.5. This is the maximum contribution for an E~ 2 signal from 
those sources not violating the diffuse limit of AMANDA. However, there 
is another production scenario which leads to a significant reduction of a 
contribution from FR-II galaxies. The synchrotron spectral index of FR-II 
ga laxies a? is correlated to th e neutrino spectral index p by p = 2 • a + 
1 (IRybicki fc Lightmanl . 1 19791 . e.g.). The mean synchrotron spectral index is 
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Fig. 11. Contribution to the diffuse neutrino background from FR-II radio 
galaxies. For an E~ 2 spectrum, it is constrained to an optical depth 
of r e ff < 0.5. A steeper spectrum E~ 2 ^ (model 2) is consistent with 
the current limit. 



a ~ 0.8 for large scale emission and therefore, the neutrino spectrum can be 
as steep as E~ 2 - 6 . This prediction is shown in Fig. 11 as model (2). In that 
case, the contribution is very low and does neither violate the current diffuse 
limit, nor is it wit hin reach of IceCube's sensitivity for t hree (or even five) 
years of operation (lAchterberg et al. (IceCube Coll.l 120071 ) . In such cL CcLSG. till 
observation of a neutrino signal from FR-II galaxies would be very difficult. If 
the spectrum were this steep for real, then all Ultra High Energy Cosmic Ray 
(UHECR) arguments would fail and so it may be more appropriate to assume 
a flatter spectrum. The radio spectral index of the hot spots are usually about 
a ~ 0.6 and so correspond to a particle spectral index of about p ~ 2.2. Using 
Fig. 11 it can be estimated that a spectrum with an index of 2.2 could already 
be explored by the three year sensitivity of IceCube. 



5.4-2 Flat Spectrum Radio Quasars 

Using a similar correlation between radio and neutrino flux as for FR-II galax- 
ies, the contribution from FSRQs to the diffuse neutrino flux can be calculated 
to 

E 2 ^- = 1.70 ■ r e// 10- 9 GeV cm" 2 s" 1 sr" 1 . (22) 
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The flux lies well beyond the current diffuse limit and even below the sensi- 
tivity of IceCube for an optical depth around r e ff = 1. Equivalent as for the 
case of FR-II radio galaxies, the upper limit to the optical depth in FSRQs 
can be determined to be 

Teff < 44 (23) 
by using the diffuse AMANDA limit. 

The effective optical depth is quite uncertain in the case of Flat Spectrum 
Radio Quasars. Just as an example, a sub-class of FSRQs are GPS, which are 
compact AGN where it is assumed that the jet runs into dense matter, yielding 
a high potential for proton-photon interaction which results in a high optical 
depth. However, other sources like TeV blazars contribute to the source class, 
being at the other end with a very low optical depth. 

With IceCube, it will be possible to examine sources with r e ff > 2.5. Con- 
sequently, neutrino emission in coincidence with the radio signal from FSRQs 
will only be observable in the near future, if the optical depth of FSRQs is 
significantly higher than the one for FR-II galaxies. 



6 Examination of source class capabilities 

In the previous sections it has been shown that the stacking method can be 
used to increase the diffuse sensitivity of high energy neutrino telescopes to 
certain source classes. 

Further catalogs published recently which are not yet part of the stacking 
search in AMANDA will be discussed in this section. The references for the 
catalogs are given throughout the text. They will be examined with respect 
to the potential neutrino signal and the effectiveness of the method. 

Also, a diffuse component from unresolved sources is expected to show devi- 
ations from a purely isotropic flux, since it is expected to follow the source 
distribution in the sky. The total signal from a certain source class to be ob- 
served by a neutrino telescope is highly dependent on the detector's field of 
view. We will generically examine source classes concerning the total flux from 
the northern respectively from the southern hemisphere. This gives a quali- 
tative examination of the capabilities of IceCube and KM3NeT which is 
planned to observe the southern hemisphere. The local supercluster for ex- 
ample is observable from the northern hemisphere and so in some cases, a 
significant fraction of the total flux in a sample comes from that hemisphere. 
It is shown, however, that there are classes with the dominant contribution in 
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the southern sky. A further constraint is the limitation of the source catalogs 
themselves: radio data are mainly given for the northern hemisphere and there 
are only a few southern identified sources in the case of very sensitive radio- 
selected samples. Here, high energy photon catalogs which are mostly provided 
by satellite experiments are much more complete in the sense of directionality. 



6.1 Additional source catalogs 



This subsection gives an overview of four additional source catalogs which 
yield high capabilities for neutrino searches with the stacking approach. The 
point source sensitivity is likely to be increased in all the cases, diffuse limits 
can be derived in two of the cases. A summary of the basic properties of these 
catalogs is given in table 3. The maximum diffusive factor for the examined 
catalogs can be determined as described in Section 4, Equ. (6), assuming that 
the total number of sources in that class can be determined. This is quite 
challenging, since a lower luminosity limit is difficult to determine. Figure 12 
shows the increase of the diffusive factor with N tot for the whole sky for all four 
examined catalogs. Table 4 reviews the parameters used for the four catalogs. 



Catalog 


energy range 


Reference 


Underlying u model 


COMPTEL 


< 100 MeV 


(KaDDadath et al.. 1996) 


( Mannheim. 1995: Stecker. 2005) 


INTEGRAL 


hard X-ray 


fBeckmann et al.. 2006) 




INTEGRAL 


soft 7- ray 


(Beckmann et al.. 2006) 




Starburst 


FIR 




fWaxman k Bahcall. 1999) 



Table 3 

Summary of source catalogs interesting to examine with respect to the neutrino 
output of the sources. 



6.1.1 Sources of Me V emission 



Neutrino flux predictions using the diffuse background as measured in hard 
X-rays by COMPTEL are based on the assumption that the cosmic diffuse 
photon flux below 100 MeV is directly correlated to the neutrino flux. About 
10 blazars could be resolved by COMPTEL underlining the assumption that 
intrinsically weaker AGN are responsible for the diffuse signal. No neutrino 
stacking analysis of the resol yed COMPTEL b lazars has bee n done y e t. Us- 
ing the same assumptions as iMannheim fll995l ) {M95-B) and ISteckerl d2005l ) 
(StSa05) would challenge such an analysis. The derivation of a diffuse limit 
from a stacking limit would however be difficult at this stage, since only a neg- 
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Catalog 


^(sources) 
south/north 


qcat 
^tot 


It/ c CJ K 




1 1 

5/6 


0.326/0.311 


01 QQ v 1 n 47 prfr/t; /Hnr 2 

0.0384/0.0199 


TNTFHR AT, 
(hard X-rays) 


-LO 

10/5 


1 7t V 10 ^ prcr /q / r*m 2 

130.08/45.1 


flQSi v If! -11 pru/s/™ 2 

0.98/4.34 


TNTROR AL 

(soft 7 rays) 


42 

23/19 


441 y ID - H prc/s/rm 2 
265/176 


85 x 10 -11 prp/s/rm 2 
1.31/0.85 


Starburst 


199 

46/153 


17000 mJy 
2260/9480 


0.906 mJy 
1.52/0.906 



Table 4 



Summary of the main parameters in the source catalogs. The total number of sources 
in a catalog is given, ^(sources). 5™/ is the integrated flux in the catalog and S wea k 
is the measured flux of the weakest source in the sample. Units are given for each 
individual source class. The second row for each class shows the same properties for 
north and south (south/north). 



* 10 4 



10 3 



10 2 



10 



1 

3 4 5 

10 10 10 

N tot 

Fig. 12. Behavior of the maximum diffusive factor with the total number of 
sources in the source class. Only a fraction of the total number of 
sources is included in the catalog, dependent on the instrument's 
sensitivity. 
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ligible fraction of the sources potentially responsible for a diffuse signal is re- 
solved into sources. Thus, the major fraction of the diffuse signal would not be 
included in a stacking analysis and the limit would not apply to diffuse predic- 
tions . The planned la unch of satellites such as GLAST (IGehrels fc Michelsonl . 
1999h and MEGA0 f lKanbach et all I2004T ) encourages to pursue the exami- 



nation of this source class, since it gives the prospect of increased statistics of 
resolved sources in the considered energy range. 



6.1.2 Soft and hard X-rays 

As mentioned before, the limit to a neutrino flux from RO SAT-detected 
sources would be improved significantly by selecting radio-weak sources in- 
stead of radio strong ones. A significant neutrino flux would not be expected 
as discussed in Section 5, but the improvement of the limit would underline 
the fact that X-ray emission from the foot of the jet is not correlated with a 
neutrino signal. 

In hard X-rays, 15 sources have been identified in the INTEGRAL0 sur- 



vey by JEM-X (IBeckmann et all 120061 ) and could be interesting to examine 



With the current analysis, HEAO-A results were used, while the more recent 
measurements by INTEGRAL are bound to be more accurate. Since a direct 
correlation between hard X-ray sources as detected by Integral and soft X-ray 
sources as seen by ROSAT does not seem to exist, an independent analysis 
would be reasonable. The emission of hard X-rays is likely to be correlated 
with TeV photon emission at least via the up-scattering of synchrotron pho- 
tons by Inverse Compton scattering. Therefore, under the assumption of a 
correlated emission of hard X-rays and TeV photons, an analysis of the hard 
X-ray signal gives an indication of the flux from high-frequency peaked BL 
Lacs. 



6.1.3 Soft 7 -rays 



A catalog of 42 AGN wit h emission in soft 7-ra ys, (20, 100) keV has been 
released by INTEGRAL ( Beckmann et all 2006 ). Here, some of the sources 
detected by OSSE^ in the energy range of (50, 150) keV could be confirmed, 
others were not seen. We suggest to use INTEGRAL data to define a source 
class for the stacking of potential neutrino sources. 



Medium Energy Gamma-ray Astronomy 

8 INTErnational Gamma-Ray Astrophysics Laboratory 

9 Oriented Scintillation Spectrometer Experiment 
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6.1.4 Starburst Galaxies 



The em ission of neutrinos from Starburst Galaxies has been suggested by lLoeb fe Waxman 



( 120061 ) (L&W in the following). It is assumed that 



(1) Relativistic protons are accelerated along with relativistic electrons, 

(2) the observed radio emission results from pion-induced electrons. The same 
pions produce neutrinos, 

(3) protons lose all their energy in pj interactions before reaching the diffu- 
sion time. 



It has been pointed out by ISteckerl (120071 ) that L&W overestimate the frac 



tion of the diffuse far infrared (FIR) flux coming from Starbursts. While it 
is ~ 23% on average, L&W assume that 100% of the detected signal comes 
from Starbursts. On the whole, the diffuse flux from Starbursts is presumable 
much lower than predicted: L&W assume that Starbursts are loss dominated, 
which means that most primaries interact and do not escape the source. This 
enhances the neutrino flux, since basically all protons lose their energy in p'-f 
or pp interactions and produce neutrinos. Observations of the spectral radio 
index of the sources (~ i^ -0,8 ) indicate however that Starbursts are in the 
diffusion limit, indicating that a negligible fraction of protons interact and 
only few neutrinos are produced. This contribution cannot be expected to be 
observed by IceCube. 

There is another possibility to expect enhanced neutrino emission from Star- 
bursts. In the past few years, it could be shown that long GRBs are typically 
connected to the explosion of Wolf-Rayet stars into a supernova Ic. These oc- 
cur preferably in St ar Forming Regions. Thus, a diffuse flux of GRBs similar 



to the prediction of IWaxman fc Bahcalll (119991 ) should originate from the di- 



rection of these galaxies. There are two different ways to normalize the diffuse 
GRB spectrum. One method is to assume that the observed keV-photon flux is 
proportional to the neutrino flux. In that case, the normalization is dependent 
on the number of observed GRBs per year. This number is strongly dependent 
on the instrument and the number is not very exact. Under the assumption 
that GRBs accelerate protons up to the highest energies, E p ~ 10 21 eV, the 
neutrino spectrum can also be normalized to the flux of ultra high energy cos- 
mic rays (UHECRs). In this case, the normalization is independent of GRB 
observations. It should be kept in mind that the spectral index of the spec- 
trum still varies from burst to burst - in the model of Waxman&Bahcall, an 
average spectral index has been used. 

It is possible to look for a neutrino signal from GRBs by stacking Star Burst 
Galaxies. This method has one advantage over a triggered-GRB search: it is 
a systematic search, since independent of GRB data. A disadvantage is that 
only nearby events are included, since the sample of Starbursts only reaches 
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out to redshifts of z — 0.07. It should however be possible to use IRASI 10 
data to identify Starburst Galaxies at higher redshift. The search for a GRB 
signal from Starbursts should be considered as a systematic search for choked 
and undetected GRBs. The sources can be selected according to their FIR- 
flux, since this is a measure of the SN rate in a Starbursts. A higher FIR flux 
indicates a high star formation rate, thus more SNe and therefore, also more 
GRBs. 



6.1.5 General approach to optimize stacking for diffuse interpretation 



In order to get diffuse limits from source stacking, it is important to choose 
source classes which have information on both resolved sources and diffuse 
background. That way, the parameters e and £ are easily and correctly de- 
termined. A good example is the EGRET catalog, where both diffuse and 
resolved emission could be proven to be correlated. Alternatively, the deter- 
mination of an upper limit to £ is possible when working with a relatively 
complete catalog of the strongest sources in the sky by using a generic num- 
ber of total sources in the source class as it is described in Section 4, Equ. (6). 
The diffuse interpretation of source classes like TeV blazars on the other hand 
is difficult, since IACTs can barely look for diffuse emission and the sensitivity 
of all-sky monitors such as MILAGROEH is not high e noug h yet to detect 
extragalactic diffuse emiss ion. Future Projects like HAWCI 12 1 ex periment as a 
successor of MILAGRO flSinnis et all [20051 : ISmith et allj2006h and CTAS 
following in the fo otsteps of H.E .S.S and MAGICI 15 1, but aiming at a 
large field of view (ITeshimal . 
of the sample. 



20061 . e.g.), will help enhancing the completeness 



In the first approach of AGN stacking in order to examine a potential excess 
in neutrinos, the number of used sources was determined by optimizing the 
significance in the detector. The optimal number of sources was in that case 
typically around ~ 10. In order to achieve the best stacking diffuse limit it 
is interesting to optimize the number of sources by taking into account £ as 
well. It needs to be tested if, by reducing £ as much as possible and taking 
the penalty of a possibly increased point source stacking limit instead, the 
diffuse limit can be increased. Also, it needs to be considered carefully if an 
optimization to a possible detection or to a limit is the most reasonable choice. 



InfraRed Astronomical Satellite 

Multiple Institution Los Alamos Gamma Ray Observatory 
High Altitude Water Cherenkov 
Cherenkov Telescope Array 
High Energy Stereoscopic System 

Major Atmospheric Gamma Imaging Cherenkov Telescope 
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6.2 Source class evolution 



For the standard analysis of muon neutrino signatures in high energy neu- 
trino detectors, the field of view is 2 it sr. The two experiments IceCube 
and KM3NeT will be observing the northern respectively the southern hemi- 
sphere in the near future and whole sky coverage is achieved by the combina- 
tion of the two experiments. 

Previously, in the stacking analysis of AMANDA data, sources with 6 < 10° 
were excluded due to the decreasing sensitivity and high muon background to- 
wards the horizon. Here, we include all sources, for both northern and southern 
hemisphere down to 5 = 0°, expecting a much better sensitivity and better 
muon rejection near the horizon for both IceCube and KM3NeT due to 
improved directional reconstruction. 

In this section, the different source classes discussed above will be examined 
with respect to their luminosity evolution. The total sample will be compared 
to the contribution from the northern and from the southern sky. For each sam- 
ple, two figures will be discussed, (a) represents the differential source counts, 
logiV(> S) versus the logarithm of the flux S. (b) shows the integral source 
evolution with the total flux of N sources, organizing the sources according to 
their strength, starting with the brightest one. Open (black) sources display 
in both cases the whole sky. Red stars consider only the southern hemisphere 
and blue triangles show the northern hemisphere contribution. 

6.2.1 Optically thick blazars - MeV to GeV emission 

Figure 13 shows the luminosity evolution of EGRET sources, for the whole 
sky as well as for the northern and southern populations. The three most lumi- 
nous sources in the EGRET sky are in the southern hemisphere. That is why 
the main contribution from GeV 7 rays is located in the southern hemisphere 
where the neutrino signal contribution is thus expected to be much higher. 
Given the restrictive limit which could already be derived from AMANDA 
data, this source sample is interesting for IceCube. Considering the source 
distribution in the sky, the class of GeV blazars is particularly useful for 
KM3NET given the high total flux. 
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F y [10 /s/cni ] N sources 

Fig. 13. (a): Number of sources in the EGRET catalog with a flux > F~. 

The complete catalog is shown as open circles. The sample has been 
divided into a sub-class of sources at the northern (triangles) and 
southern (stars) hemisphere, (b): Total luminosities of all sources 
< N, starting to sum up with the strongest source and successively 
adding the next luminous source. 

The COMPTEL Catalog of sources with E < 100 MeV is displayed in Fig. 14. 
In this case, the main contribution lies in the northern hemisphere. The small 
number of only 11 identified COMPTEL sources make a diffuse interpretation 
of a possible limit difficult. Satellites like MEGA would make the investiga- 
tion even more interesting. The 11 identified sources are however still useful to 
investigate with respect to a point source signal. The maximum diffusive factor 
lies around £ max ~ 300 for 10,000 sources and £ max ~ 3000 for 100,000 sources 
in the class. Depending on the steepness of the evolution function, this factor 
can be significantly smaller, so that doubling the number of sources could al- 
ready help to draw a conclusion about a diffuse limit. 



northern sources 
southern sources 



6.2.2 Soft ^-rays 



A cat alog of soft 7-ray-detected AGN is given by INTEGRAL (IBeckmann et al. 



20061 ) . 19 of 42 sources are in the northern sky, among the 23 southern sources 
are the three most luminous ones. The catalog is a good candidate for dif- 
fuse interpretation of the neutrino results, given the relatively high number 
of sources. The luminosity evolution is displayed in Fig. 15. It can be seen 
that the evolution is still rising and that there is still a significant fraction 
of signal missing. This is mirrored in the numerical value of the maximum 
diffusive factor which is calculated to £ max ~ 20 for N tot = 10, 000 sources 
and £ max ~ 200 for Ntot — 100, 000 sources. These relatively large numbers 
increase a potential limit on the neutrino flux from INTEGRAL sources by 



about one order of magnitude 16 



16 In Section 4 it is already discussed that even with a source class of 10 5 sources, 
the strongest 10 4 AGN make up the dominant contribution. Thus £ max = 20 is the 
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Fig. 14. (a): Number of sources in the COMPTEL catalog with an energy 
density > p e . (b): Total energy density of all sources < N, starting 
to sum up with the strongest source and successively adding the next 
luminous source. 
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Fig. 15. (a): log7V(> 5) - logS-plot for INTEGRAL sources - total cat- 
alog (circles), only northern sources (triangles) and only south- 
ern sources(stars). (b): Integral flux versus number of contributing 
sources. 

6.2.3 Soft and hard X-rays 

Figure 16 shows the luminosity evolution for RO SAT-detected sources. 84 
sources are in the sample totally, of which 35 are northern and 49 are south- 
ern sources. The three most luminous sources are in the northern hemisphere. 



more realistic value in this case. 
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Fig. 16. (a): Number of sources in the ROSAT catalog with an energy density 
> p e . (b): Total energy density of all sources < N, starting to sum up 
with the strongest source and successively adding the next luminous 
source. 



However, these are objects which have such a large flux, since they are ex- 
tremely close to Earth and not because of their high intrinsi c luminosity. This 



i s why these sources were excluded in the stacking analysis of lAchterberg et al. (IceCube Coll. fc P. L. I 

( 120061 ) . Without these sources, the contributions from both hemispheres are 

comparable. 
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Fig. 17. (a): logiV(> S) - log S-plot for INTEGRAL sources in hard X-rays 
(2, 10) keV - total catalog (circles), only northern sources (triangles) 
and only southern sources (stars), (b) total flux of N contributing 
sources, starting with the brightest sources and going to weaker fluxes 
with higher N. 
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A sample of hard X-ray AGN is available from the INTEGRAL satellite, 
see Fig. 17. INTEGRAL detected 15 sources in the energy band of 2 keV 
to 10 keV. Among the 10 sources in the southern hemisphere are the two 
strongest sources in the sample. There are only 5 northern sources and the 
main contribution comes from the southern hemisphere. The INTEGRAL 
sample can improve the stacking analysis of hard X-ray sources which has been 
done with HEAO-A information. Only three northern sources were reported 
from HEAO-A. For a diffuse interpretation, there are still too few sources in 
the sample, though. 



6.2.4 FR-I/FR-II 

The catalog FR-I and FR-II sources is restricted to values of 5 > —10°. There- 
fore, only few southern sources (7 FR-I and 14 FR-II galaxies) are in the com- 
plete sample which is seen in Fig. 18 for FR-I galaxies and in Fig. 19 for FR-II 
galaxies. The flux is totally dominated by the northern hemisphere. 
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Fig. 18. (a): Number of sources > S for FR-I galaxies, (b) Total flux versus 
number of sources N. 



6.2.5 CSS /GPS 

The problem for the catalog of CSS and GPS sources is similar to the situation 
of the FR-I and FR-II catalog. CSS have been selected at 5 > 10°, while GPS 
include data with S > —25°. Therefore, there are no southern sources in the 
case of CSS. Figure 20 shows the GPS sample. 8 southern and 20 northern 
sources have been identified with the main contribution from the northern 
hemisphere. 
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Fig. 19. (a): Number of sources > S for FR-II galaxies, (b) Total flux versus 
number of sources N. 
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Fig. 20. Giga-Hertz-Peaked Sources (GPS) (|0 'Pea fc Bauml . fl997l l . (a) is the 
differential flux evolution, (b) shows the total flux versus number of 
sources N. 



6.2.6 QSOs 



In the sample of QSOs as presented by ISanders et al.l (119891 ). there are only 
3 sources in the southern hemisphere with a measured flux at the selection 
wavelength, 60 /zm. The sample selection of the sample was done for 8 > —15°, 
which only leaves a small window on the southern sky. While this source class 
is well-suited for the analysis of the northern hemisphere, there is too little 
data available in the southern hemisphere. 
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6.2.7 Starburst Galaxies 



The catalog of Starburst Galaxies includes 199 sources. While the most lu- 
minous sources are in the southern, the total flux is higher in the northern 
hemisphere, since 153 of 199 sources are located north. A stacking analysis 
has good potential for both hemispheres. 




Fig. 21. (a): Number of sources in the Starburst catalog with a flux > F~. (b): 
Total flux of all sources < N, starting to sum up with the strongest 
source and successively adding the next luminous source. 



7 Conclusions and Outlook 



In this paper, limits to the diffuse flux from seven different AGN source classes 
could be derived from point source stacking limits from AMANDA. There is 
in many cases an additional constraint on flux models at the lowest energies: 
while the general diffuse limit is restricted to the energy band of 10 4,2 GeV 
to 10 6 ' 4 GeV, the stacking limits reach down to 10 3 GeV. In many cases, the 
diffusive factor can conservatively only be estimated. It is expected, that the 
additional diffuse component is much smaller than the values derived here as 
upper limits. A reduction of the diffusive factor would result in the improve- 
ment of the stacking diffuse limits. 



Neutrino flux predictions normalized to the diffuse EGRET signal above 
100 MeV use the as sumption that the diffuse componen t is pr o duced by un 



resolv ed AGN, see ( iMannheim et all l200ll ; iMannheiml . Il995l ; iMiicke et al. 
20031 ). Thus, the EGRET stacking limit applies to such calculation s . The 



limit falls short of the outline of two of the models, (IMannheim et all 12001 
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Mannheiml . Il995l ) . The model of proton-proton interactions normalized to the 
EGRET diffuse results violates the limit at low energies around 10 3 GeV. At 
the highest energies, the stacking diffuse limit is more restrictive than the gen- 
eral diffuse limit and reaches below sensitivities of the maximum contribution 
from proton-photon interactions. A detection of a flux from optically thick 
blazars should be visible with second generation neutrino telescopes such as 
IceCube. 

The neutrino flux models from X-ray detected AGN and from FR-II galaxies 
can be constrained by the general diffuse limit. 

In the case of X-ray AGN, the limit strongly disfavors a hadronic model. It is 
an order of magnitude below the flux predictions. According to the strongly 
restrictive limit, neutrinos are not produced at the foot of AGN jets in coin- 
cidence with X-rays. Thus, the diffuse emission as measured by ROSAT is 
likely to be due to Inverse Compton radiation. 

Analyzing the correlation between the radio and neutrino emission from FR- 
II galaxies and blazars, the limit was used to constrain the optical depth of 
the sources. For FR-II galaxies, r e // < 0.5 could be derived, while the upper 
limit for blazars is given as r e // < 44. Thus, the detection potential for such 
a source class in IceCube is very high, since within three years, IceCube is 
sensitive to sources of optical depth with r e // > 0.024 (FR-II) resp. r e // > 2.1 
(blazars) . 

It was shown that a general diffuse analysis gives a high discovery potential 
for hidden TeV sources: While TeV photons are absorbed at high redshifts, 
neutrinos propagate freely. An upper limit to the neutrino flux from photon- 
resolved TeV blazars is determined to be 

$ DL resoivedTeV = 1-37 ■ I0~ 9 GeV cm" 2 s" 1 sr" 1 . (24) 



The considerations above show that an investigation of a neutrino signal can 
be used to help determining different intrinsic parameters of the source type 
considered. It could be shown that a stacking analysis yields valuable informa- 
tion on the diffuse contribution from theses sources. The detection probability 
is enlarged significantly in such an approach and there is a high detection 
potential with IceCube and KM3NeT. A cascade stacking analysis by these 
neutrino detection arrays should be considered in order to incre ase the sensi- 



tivity to ultra high energy neutrino fluxes as they are described by lMiicke et al 



fl2003|). 



Apart from the catalogs which have already been used in the stacking analy- 
sis of AMANDA, further source classes are investigated here. The possibility 
of a stacking analysis of COMPTEL sources is discussed in this paper as 
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a first approach to examine optically thick sources with photon emission of 
E < 100 MeV. Although no stacking diffuse limit can be derived yet due to the 
small percentage of resolved sourc es, future experiments such a s GL AST at an 



energy range o f (0.02,300) GeV flGehrels fc Michelsonl . I1999T ) and MEG A at 



(0.4, 50) MeV (IKanbach et al.l . 120041 ) give hope to resolve many more sources 
in the MeV range. Another interesting source class is the catalog of Starburst 
Galaxies. With the correlation between starforming regions and long duration 
gamma ray bursts, an enhanced neutrino signal from Starbursts can be ex- 
pected. The prospects of IceCube and KM3NeT are different considering 
source stacking. Most catalogs show a bias to one of the hemispheres. Both 
neutrino detection arrays can be used to extract complementary information 
about the different source classes. 

The upper limits derived in this paper constraint several prevailing neutrino 
flux models. This underlines the necessity of calculation on the basis of new de- 
velopments within astroparticle physics and numerical approaches. A unified 
model concerning acceleration processes in AGN would give the opportunity 
to examine a potential neutrino signal with respect to different AGN classes 
within the same framework. The numerical results could be applied to indi- 
vidual sources and to diffuse flux measurements. The interaction between the 
observation of resolved sources and diffuse photon and proton components of 
CRs is of high significance with respect as it could be shown in this paper. 
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